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PURPOSE:  This  Coastal  and  Hydraulics  Engineering  technical  note  (CHETN)  describes  an 
automated  method  to  develop  a  synthetic  unit  hydrograph  of  a  small  watershed  using  readily 
available  spatial  data  products,  namely  Digital  Elevation  Model  (DEM),  and  land  cover  data.  A 
Python  script  is  utilized  within  ArcGIS  (Environmental  Systems  Research  Institute  2011)  to 
quickly  and  efficiently  determine  how  the  spatial  features  of  the  watershed  affect  the  runoff  of 
the  basin  and  therefore  the  unit  hydrograph  at  the  outlet  of  the  basin. 

BACKGROUND:  Rainfall-runoff  response  within  a  watershed  is  a  core  consideration  of 
hydrologists.  The  use  of  unit  hydrographs  as  a  way  to  analyze  the  rainfall-runoff  response  was 
first  proposed  by  Sherman  (1932).  Since  then,  many  interpretations  of  the  unit  hydrograph  have 
been  developed  that  assume  runoff  is  routed  through  a  series  of  linear  reservoirs,  including  the 
following:  instantaneous  unit  hydrograph  method  (Nash  1957),  geomorphic  instantaneous  unit 
hydrograph  method  (Rodriguez-Iturbe  and  Valdes  1979),  and  the  Clark  unit  hydrograph  (CUH) 
method  (Clark  1945).  The  effectiveness  of  the  unit  hydrograph  is  well  evident  in  hydrologic 
practice,  with  the  popularity  of  the  CUH  method  having  been  increased  particularly  due  to  its 
inclusion  in  the  U.S.  Army  Corps  of  Engineers  (USACE)  Hydrologic  Engineering  Center  (HEC) 
suite  of  watershed  models,  including  HEC-1  (Eeldman  1995)  and  the  Hydrologic  Modeling 
System  (Scharffenberg  and  Eleming  2006). 

Within  gauged  basins,  the  unit  hydrograph  can  be  developed  using  observed  storm  events 
(precipitation  and  runoff  data  required).  In  basins  with  little-to-no  observed  data,  a  synthetic  unit 
hydrograph  can  be  developed.  The  CUH  is  a  synthetic  unit  hydrograph  that  is  developed  using  a 
time-area  method.  To  develop  a  hydrograph  using  the  CUH  method,  three  parameters  are 
required:  time  of  concentration  {Td),  storage  attenuation  coefficient  (R),  and  a  time-area 
histogram  (TAH)  of  the  basin.  Tc  is  a  measure  of  the  time  that  it  takes  for  water  to  flow  from  the 
most  remote  part  of  the  basin  to  the  outlet.  Theoretically,  from  the  time  rainfall  begins,  all  of  the 
runoff  from  the  entire  basin  should  be  contributing  to  the  outlet  when  the  time  reaches  Tc.  R 
represents  the  effects  of  storage  on  the  timing  of  the  outflow  from  the  basin.  TAH  is  a  means  to 
represent  the  area  within  the  basin  that  drains  to  the  outlet  at  given  time  intervals  and  is  the  focus 
of  this  CHETN. 

Previous  methods  have  used  GIS  to  help  create  a  CUH  for  a  basin  (Eabadie  2014;  Maidment  et 
al.  1996;  Maidment  1993).  To  help  improve  implementation  times  during  emergency  scenarios, 
this  document  shows  how  a  TAH  can  be  automatically  derived  using  spatially  heterogeneous 
parameters  based  on  topographic  (derived  from  DEM)  and  land  cover  characteristics.  Because 
the  TAH  is  calculated  based  on  routing  times,  Tc  is  intrinsically  included  within  TAH.  For  the 
purposes  of  this  document,  any  delay  in  outflow  due  to  R  is  also  assumed  to  be  included  within 
the  flow  routing  scheme.  A  simple  Python  script  was  created  that  automates  the  generation  of  the 
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CUH.  The  Python  script  (Appendix  A)  runs  within  ArcGIS  and  generates  the  TAH  by  calling 
ArcGIS  functions  and  a  simple  executable,  whose  code  is  shown  in  Appendix  B.  The  format  of 
this  CHETN  is  to  show  the  methodology  on  how  the  Python  script  works  in  a  step-by-step 
example.  The  purpose  of  this  CHETN  is  to  show  how  readily  available  spatial  data  can  be  used 
within  GIS  to  create  a  simple,  yet  accurate,  spatially  discretized  CUH. 

METHODS/EXAMPLE:  The  test  location  selected  for  this  CHETN  is  the  New  River  basin  in 
coastal  North  Carolina  (Eigure  1).  The  basin  is  approximately  211  km  and  ranges  in  elevation 
between  0.55  m  and  36.80  m.  As  presented  in  Table  1,  the  area  is  predominately  evergreen 
forest,  cultivated  crops,  and  woody  wetlands.  This  test  location  was  selected  because  it  is  near 
Marine  Corps  Base  Camp  LeJeune,  the  site  of  a  US  ACE  test  project  aimed  at  deployment  of 
fast-response  hydrologic  models. 
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Table  1.  NLCD  2006  identification  number,  classification,  associated 
overland  flow  conveyance  (K)  values  adapted  from  McCuen  (1998),  and 
runoff  coefficient  (Cr)  values  taken  from  the  American  Society  of  Civil 
Engineers  (ASCE)  and  Water  Pollution  Control  Federation  (WPCF)  (1982). 


NLCD  2006 

ID 

Land  Cover 

%  Area 

K 

C,  (fraction) 

11 

Open  Water 

0.07 

25.00 

1.00 

21 

Developed,  Open  Space 

3.84 

5.10 

0.20 

22 

Developed,  Low  Intensity 

1.02 

10.20 

0.30 

23 

Developed,  Medium  Intensity 

0.25 

15.30 

0.40 

24 

Developed,  High  Intensity 

0.01 

20.40 

0.50 

31 

Barren  Land 

0.39 

2.00 

0.15 

41 

Deciduous  Forest 

0.11 

1.60 

0.15 

42 

Evergreen  Forest 

20.84 

1.60 

0.15 

43 

Mixed  Forest 

3.96 

1.60 

0.15 

52 

Shrub/Scrub 

8.90 

7.00 

0.15 

71 

Herbaceuous 

2.86 

7.00 

0.15 

81 

Hay/Pasture 

0.10 

3.50 

0.15 

82 

Cultivated  Crops 

31.82 

2.00 

0.15 

90 

Woody  Wetlands 

25.00 

1.60 

0.15 

95 

Emergent  Herbaceuous  Wetlands 

0.81 

1.60 

0.15 

The  following  is  a  step-by-step  example  on  how  a  CUH  can  be  easily  derived  at  a  point  along  the 
river  within  the  New  River  basin.  The  data  used  within  the  example  are  freely  and  readily 
available  nationwide,  so  the  processes  and  methods  shown  can  be  implemented  easily  at  another 
location.  An  example  of  how  flow  can  be  calculated  using  the  CUH  method  is  also  shown  and 
compared  to  observed  data. 

Step  1.  Obtaining  spatial  data.  The  1/3  arc  second  (~9.645  m)  National  Elevation  Dataset 
(Gesch  et  al.  2002)  over  the  project  area  was  obtained  from  the  United  States  Geological  Survey 
(USGS)  (http://nationalmap.sovA.  A  land  cover  map  from  the  National  Land  Cover  Database 
2006  (NLCD  2006)  (Lry  et  al.  2011)  was  obtained  from  the  same  website.  Both  datasets  were 
reprojected  onto  a  UTM  Zone  18N  projection  (the  local  UTM  projection). 

Step  2.  Create  flow  direction  and  fiow  accumuiation  rasters;  specify  outiet.  Llow 
direction  (LDIR)  and  flow  accumulation  (LAC)  rasters  were  created  over  the  entire  area  using 
the  “Hydrology  Toolbox”  within  ArcGIS.  Lor  each  cell  in  the  domain,  LDIR  specifies  the 
direction  of  flow  from  a  cell  to  its  neighboring  cell,  and  LAC  indicates  how  many  up-slope  cells 
hydrologically  contribute  to  the  cell.  Based  on  the  LAC  raster,  a  point  shapefile  was  created  at 
the  outlet  of  the  basin  being  examined.  Ligure  2  shows  the  LAC  raster  as  well  as  the  location  of 
the  outlet  selected.  Steps  1  and  2  are  typical  processes  carried  out  prior  to  any  hydrologic 
investigation  and  therefore  were  not  included  in  the  Python  script.  The  remaining  steps  are 
carried  out  within  the  Python  script  shown  in  Appendix  A. 
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Step  3.  Creation  of  TAH.  There  are  two  products  within  this  step:  (1)  a  raster  that  represents  the 
time  it  takes  each  cell  within  the  basin  to  drain  to  the  outlet  and  (2)  a  TAH  derived  from  the  raster. 
To  complete  this,  several  subtasks  must  be  completed.  These  tasks  are  automated  within  a  Python 
script  shown  in  Appendix  A  and  a  C++  code  shown  in  Appendix  B.  Tasks  1-7  in  the  following  list 
provide  a  description  of  what  is  simulated  within  the  Python  script,  while  Tasks  8  and  9  provide  a 
brief  description  of  how  the  two  products  are  ultimately  created  using  a  C++  code. 

Task  1  (Lines31-41):  The  watershed  is  delineated,  and  a  watershed  boundary  is  created  (Figure  2). 
The  elevation  and  flow  accumulation  rasters  are  clipped  to  the  watershed  boundary.  A  stream  grid 
is  created  by  assuming  that  any  FAC  cell  that  drains  an  area  greater  than  5  km  (cell  value  of 
53,744  at  a  9.645  m  resolution)  is  considered  a  stream. 

Task  2  (Lines  44-48):  The  slope  of  each  cell  is  calculated  using  the  “Slope”  function  within 
ArcGIS.  The  slopes  are  calculated  for  an  area  larger  than  the  watershed  boundary  and  then  are 
clipped  to  the  watershed  boundary.  The  clipped  slope  raster  is  initially  calculated  as  a  percentage 
with  all  slopes  less  than  0.001%  being  assigned  a  value  of  0.001%.  The  slope  raster  is  then 
converted  into  a  decimal  format  raster  (OutSlope  Dec).  A  third  raster  is  also  created  to  represent 
the  slope  (decimal  format)  in  only  the  channel  cells  (Slope  Chan). 

Task  3  (Lines  51-58):  The  land  cover  raster  is  clipped  to  the  watershed  boundary  and  is  set  to 
match  the  spatial  extent  and  resolution  of  the  DEM  raster.  Flow  conveyance  values  {K)  are  then 
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mapped  to  a  new  raster  (kfactor)  based  on  the  land  eover  elassifieation  and  the  assoeiated  K 
value.  Table  1  lists  the  values  of  K  used  in  this  projeet,  whieh  are  based  on  MeCuen  (1998)  and 
Nieklow  et  al.  (2006). 


Task  4  (Lines  61-63):  The  time  it  takes  for  water  to  flow  overland  within  eaeh  eell 
(timeoverland)  is  ealeulated  based  on  an  approximation  of  the  Manning’s  equation: 


timeoverland 


1mm 

1 

60  sec 

0.3048  m 

S-0.5 


0.05A7 [C  Ky^  S-°-^ 


where  K  is  taken  from  kfactor,  S  is  the  slope  of  the  cell  (decimal,  taken  from  OutSlope  Dec),  and 
timeoverland  is  in  min-m'\  Values  of  K  are  approximations  that  incorporate  channel/flow 
geometry  and  roughness.  It  is  also  realistic  to  assume  that  K  varies  based  on  scale.  A  calibration 
parameter  (Q  is  included  to  manipulate  the  value  of  K  for  a  given  basin  and  is  set  to  0.3  for  this 
example,  which  effectively  increases  the  time  it  takes  for  water  to  flow  over  each  cell. 

Task  5  (Lines  66-69):  A  raster  representing  the  time  of  travel  within  each  channel  cell 
(timechannel,  min  m'^)  is  calculated  using  the  same  approach  as  for  overland  flow,  with  the 
exceptions  that  S  is  taken  from  Slope_Chan,  and  a  static  value  of  25.0  is  used  for  K.  In  reality, 
the  K  value  likely  varies  with  each  segment  of  stream,  but  for  this  simple  example,  it  is  left  as  a 
constant  for  all  channel  cells.  Similar  to  Task  4,  a  channel  calibration  parameter  is  also  employed 
to  manipulate  the  value  of  K  and  is  also  set  to  0.3  for  this  example. 

Task  6  (Lines  72-74):  A  raster  called  traveltime  is  created  that  represents  the  travel  times 
(min-m’')  of  both  the  overland  flow  and  channel  flow.  If  a  cell  is  classified  as  a  channel,  the 
value  from  timechannel  is  used,  otherwise  the  value  is  taken  from  timeoverland.  After 
calculation,  the  traveltime  raster  is  converted  into  an  ASCII  format  {traveltime. asc)  for  use  in 
Tasks  7-9. 

Task  7  (Lines  77-78):  Although  the  time-per-distance  for  water  to  travel  over  each  cell  has  been 
determined  in  Task  6,  the  path  that  water  takes  typically  depends  on  topography  and  can  be 
approximated  as  the  path  it  would  follow  along  increasing  cell  values  of  the  FAC.  The 
accumulated  time  for  water  to  travel  a  hydraulically  driven  path  from  each  cell  to  the  outlet  is 
required  to  create  a  TAH.  A  function  within  ArcGIS  that  can  complete  this  task  does  not  exist. 
Therefore,  a  program  was  written  in  C++  to  complete  this  task  and  is  shown  in  Appendix  B. 

Task  8  (Appendix  B):  A  simple  executable  is  called  from  the  Python  script.  The  executable  was 
developed  to  determine  the  time  it  takes  for  each  cell  to  flow  to  the  outlet  based  on  the  flow 
pathway  exhibited  by  the  FAC  raster  and  the  time  to  flow  across  a  given  cell  in  Task  6.  The 
program  routes  the  path  water  takes  from  each  cell  within  the  watershed  to  the  outlet  by 
following  a  path  of  increasing  FAC  values.  The  time  (min)  it  takes  to  travel  across  a  single  cell  is 
simply  the  distance  across  the  cell  (m)  multiplied  by  the  value  within  the  traveltime. asc  raster 
(min-m’').  The  total  time  (min)  it  takes  each  cell  to  route  to  the  outlet  is  the  summation  of  the 
times  along  the  pathway.  This  summation  of  the  time  is  recorded  in  the  totaltime.asc  raster  and  is 
shown  for  the  New  River  basin  in  Figure  3.  To  clarify,  the  totaltime.asc  raster  specifies  the  time, 
in  minutes,  it  takes  for  water  to  travel  from  that  cell  to  the  outlet  of  the  basin. 
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N 


I  I  4,320  -  5.040 

I  I  5,040  -  5,760 

I  I  5.760  -  6,480 

hi  t  6,480  -  7,200 

7,200  -  7,920 
7,920  -  8,640 
8.640  + 
if  Outlet 


2.5  5  Kilometers 


Figure  3.  Time  in  minutes  for  each  ceii  to  drain  to  outiet  {totaltime. asc). 


Task  9  (Appendix  B):  A  histogram  based  on  the  totaltime. asc  raster  shown  in  Task  8  can  be 
generated  using  a  simple  plotting  program,  such  as  Excel.  Figure  4  shows  such  a  histogram  using 
6  hr  intervals,  which  is  the  6  hr  TAH. 

Step  4.  Flow  Simulation.  Two  methods  are  available  to  determine  flow  based  on  a  given 
precipitation.  For  uniform  precipitation,  a  unit  hydrograph  is  easily  generated  using  the  S-curve 
method  and  the  TAH,  where  a  standard  unit  hydrograph  can  be  derived  for  a  given  time  interval 
of  rainfall  and  used  with  the  rainfall  data,  or  a  normalized  unit  hydrograph  can  be  derived  for  use 
in  conjunction  with  peak  flow  calculations  (Bras  1990;  Brutsaert  2005).  Figure  5  shows  an  S- 
curve  with  a  normalized  unit  hydrograph.  The  second  method,  which  is  used  in  this  document, 
uses  spatial  precipitation  data  in  conjunction  with  the  totaltime. asc  raster  to  determine  how  long 
the  rain  in  each  cell  takes  to  get  to  the  outlet.  This  second  method  takes  a  Fagrangian  approach  to 
simulate  flow  because  it  follows  rainfall  through  each  cell  to  the  outlet.  Both  methods  described 
assume  the  principle  of  superposition  to  simulate  outflow. 
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TESTING:  Two  peak  flows  from  the  New  River  near  Gum  Branch  streamgage  occurred 
between  15  August  and  15  September  2011  when  multiple  large  rain  events  occurred  over  the 
study  area.  This  time  period  will  serve  to  show  how  the  CUH  method  described  can  be  used  to 
simulate  flow. 

Next-Generation  Radar  (NEXRAD)  precipitation  data  in  shapefile  format  were  obtained  from 
the  National  Weather  Service  (htW://www.srh.noaa.sov/ridse2/RFC  Preciy/).  The  data  were 
reprojected,  clipped  to  the  area  of  interest,  and  converted  into  raster  format  using  the  inverse 
distance  weighted  (IDW)  function  within  ArcGIS.  Flow  data  from  the  New  River  near  Gum 
Branch  streamgage  were  obtained  from  the  USGS  {http: //water data. usss. sov/nwis/rt). 

Using  a  simple  programming  script  (not  shown),  the  effective  precipitation  (Pep)  was  applied  to 
the  totaltime. asc  raster  to  generate  a  hydraulic  response  from  each  hour  of  precipitation.  By 
combining  these  hourly  responses  through  the  assumption  of  superposition,  a  simulated  outflow 
hydrograph  is  created.  Pgff  represents  the  fraction  of  precipitation  that  contributes  to  immediate 
runoff,  and  is  calculated  as 


Peff=PC. 


where  P  is  the  precipitation  over  an  hour  time  period  for  a  given  cell  from  the  NEXRAD  data 
and  Cr  is  the  runoff  coefficient  (Table  1).  Values  of  C,-  vary  widely  and  are  typically  calibrated 
parameters.  The  values  of  C,.  used  in  this  project  are  based  on  land  cover  classifications  and 
ASCE  and  WPCF  (1982).  Developed  areas  are  based  on  residential  and  neighborhood  values, 
while  all  other  areas  use  a  value  of  0.15,  which  is  the  approximate  average  of  the  values  for 
lawns,  unimproved  areas,  and  parks. 

The  time-series  of  hourly  precipitation  (volume,  m  )  applied  to  the  watershed  is  shown  with  the 
simulated  (Qsim)  and  observed  flow  {Qohs)  in  Figure  6.  Qsim  for  the  first  event  (26  August-4 
September)  peaked  approximately  7  hr  before  Qobs  and  overestimated  the  peak  by  approximately 
5.51  m^  s'^  Qsim  peaked  approximately  9  hr  after  Qobs  during  the  second  event  (7-10  September) 
and  underestimated  the  peak  by  approximately  4.41  m^  s'^  The  duration  of  both  events  matched 
reasonably  well  when  compared  to  Qobs-  This  highlights  potential  areas  of  improvement  in  both 
the  inclusion  of  baseflow  estimation  and  the  methods  used  to  calculate  the  time  to  flow  across 
each  cell.  For  both  events.  Table  2  shows  the  simulated  and  observed  peak  flows,  mean  bias 
error  (MBE),  root  mean  squared  error  (RMSE),  and  Pearson’s  correlation  coefficient  (R). 


MBE 


n^^J2iQsim-Qobs) 

1=1 


RMSE 


0.5 


Table  2.  Flow  event  statistics  (aii  in  units  of  s'^  except  R, 

which  is  0-1). 

Peak  Qobs 

Peak  Qsio, 

MBE 

RMSE 

R 

Flow  Event  1 

30.30 

35.81 

-1.73 

3.77 

0.95 

Flow  Event  2 

10.34 

5.93 

-1.14 

2.19 

0.84 
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Figure  6.  Test  events:  Precipitation  (soiid  biack  iine,  top),  Observed  Fiow  (soiid  biack  iine,  bottom),  and 
Simuiated  Fiow  (dotted  biack  iine). 


CONCLUSIONS:  This  CHETN  describes  in  a  step-by-step  procedure  how  to  automatically 
develop  a  time-area  unit  hydrograph  based  only  on  elevation  and  land  cover  data.  The  CHETN 
also  provides  the  programming  code,  which  can  be  used  by  anyone  with  low-to-intermediate 
ArcGIS  and  programming  skills.  It  was  also  shown  that  the  time-area  unit  hydrograph  methods 
can  in  turn  be  used  to  accurately  simulate  outflow  from  a  small  basin  when  precipitation  data  is 
known.  Through  an  example,  distributed  rainfall  data  was  applied  to  the  generated  synthetic  unit 
hydrograph  to  produce  accurate  results.  Although  the  methods  presented  are  simple,  further 
refinement  can  be  introduced  into  the  framework  to  improve  accuracy  in  time-to-drain  estimates 
and  effective  precipitation  estimation.  Although  this  project  was  developed  to  assist  in  a  specific 
military  need,  the  methods  may  be  applicable  to  other  modeling  frameworks  that  use  a  unit 
hydrograph  method. 

POINTS  OF  CONTACT:  This  CHETN  was  prepared  by  the  Military  Engineering  Program 
through  the  Austere  Entry  Assessment  project.  The  POC  for  technical  inquiries  is  Mike  Eollum 
(Michael. L.Follum(a),usace. army. mil).  This  technical  note  should  be  referenced  as  follows: 

M.  Eollum.  2015.  Automated  method  to  develop  a  Clark  synthetic  unit 
hydrograph  within  ArcGIS.  ERDC/CHL  CHETN-IV-104.  Vicksburg,  MS:  U.S. 

Army  Engineer  Research  and  Development  Center.  An  electronic  copy  of  this 
CHETN  is  available  from  http: //chi. erdc.  usace. army. mil/chetn. 
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Appendix  A:  Python  Program  Code  Executed  through  ArcGIS 


1  t  In^ort  arcpv  module 

2  import  arcpy 

3  from  arcpy  import  env 
^  from  arcpy. sa  Uigport  • 

5  import  arcpy. cartography  as  CA 

6  import  arcpy. management  as  W 
Import  arcpy. analysis  as  ANL 

8  arcpy  .CheclcOutExtension  ("spatial") 

9  arcpy .env. workspace  »  "CUH_Workspace . gdb" 

10  arcpy .env. overwriteOutput  *  True 

11 

12  MF  «  "C:\\Work_Folder" 

13  Outlet  ■  "C ; /Work_Folder/NR_GumBranch_Gage . shp" 

14  Big_DEM  -  "C:/Work_Folder/CL_DEM_10m_cl.img" 

15  Big_FAC  -  "C:/Work_Folder/CL_FAC.img" 

16  Big^FDIR  »  "C:/Work_Folder/CL_FDIR.  img" 

17  NLCD  -  "C:/Work_Folder/NC_NLCp_2006_clipped.tif" 

IS 

19  #  Local  variables: 

20  Matershed_Poly  ■  MF  +  "Watershed_Poly .shp" 

21  FAC  »  MF  +  "\\"  4-  "FAC" 

22  FAC_ASC  »  MF  +  "W"  +  "FAC_ASC.asc" 

23  kfactor  “  MF  +  "\\"  +  "kfactor" 

24  timeoverland  *  MF  +  "\\"  +  "timeoverland" 

25  timechannel  «  MF  4-  "\\"  +  "timechannel" 

26  traveltime  *  MF  +  "\\"  +  "traveltime" 

27  traveltimeasc  *  MF  +  "\\"  +  "traveltime.asc" 

28  LC_Haskasc  «  MF  +  "\\"  +  "LC_Mask. asc" 

29 

30  #  Delineate  a  Watershed  and  Create  one  with  a  mask  watershed  (Watershed_l) 

31  outSnapPour  »  SnapPourPoint (Outlet,  Big_FAC,  20,  "FID") 

32  outWatershed  *“  Watershed (Big_FDIR,  outSnapPour) 

33  Water5hed_l  *  Con (outWatershed  >»  0,  1) 

34  arcpy  .RasterToPolygon___conver Sion  (Watershed^!,  Watershed_Poiy,  "NO^SIMPLIFY",  "VALUE") 

35 

3€  Extract_DEM  **  ExtractByMask (Big_DEM,  Watershed_Poly) 

37  Extract_FAC  *  ExtractByMask (Big_FAC,  Watershed_Poly) 

38  Extract^FAC . save (FAC) 

39  arcpy .RasterToASCII_conversion(Extract_FAC,  FAC^ASC) 

40  StrmVal  -  53744.0  *53744  is  the  number  of  cells  within  5  km2 

41  ChanGrid  -  Con (Extract_FAC  >  StrmVal,  1) 

42 

43  *  Calculate  Slope  of  the  Basin 

44  OutSlope  =  Slope (Big_DEH,  "PERCENT_RISE",  1) 

45  OutSlope_Clip  *»  ExtractByMask (  OutSlope,  Watershed_Poly  ) 

46  OutSlope_Percent  ■  Con (OutSlope^Clip  <*  0.001,  0.001,  OutSlope) 

47  OutSlope_Dec  *  Divide (OutSlope_Percent,  100.0) 

4S  Slope  ^an  *  Times  (OutSlope  D^,  ChanGrid) 

49 

50  *  Land  Use  Grid 

51  tempEnvironmentO  »  arcpy. env. cellSize 

52  arcpy . env. cellSize  *  "MINOF" 

53  L0__0  “  ExtractByHask  (NLCD,  Watershed^!) 

54  arcpy. env. cellSize  »  ten^EnvironmentO 

55  LU_1  “  Con(LU_0  “  11,25.0,  Con(LU_0  =* 

56  Con(Ln_0  =  41,1.6,  Con(LU_0  = 

57  Con(LU_0  =  81,3.5,  Con(LU_0  = 

58  LU_1 . save (kfactor) 

59 

60  *  Calculate  Time  (min)  and  Velocity  (m/s)  of  Overland  Flow 

61  OF^Calib  •  0.3  #Overland  Flow  Calibration  Factor 

62  OV_2  “  0.0547  "  (1 . 0/ (OF_Calib*Ln_l) )  •  Power (OutSlope_Dec,  -0.5) 

63  OV_2 . save ( timeoverland) 

64 

65  *  Get  the  Slope  of  the  Channel  and  a  Time  of  Travel  within  the  Channel  (minutes) 

66  K_Channel  =  25.0  tConveyance  of  channel,  assumed  to  be  25 

67  CH_Calib  =0.3  IChannel  Flow  Calibration  Factor 

68  Chanel  *  0.0547  *  (1/ (K_Channel*CH_Calib) )  *  Power (Slope_Chan,  -0.5) 

69  Chanel. save (timechannel) 

70 

71  #  Compute  Travel  Time 

72  Time_0  *  Con (IsNull (Chan_l) , 0V_2 , Chan_l) 

73  Time^O.save (traveltime) 

74  arcpy .RasterToASCII^conversion (Time^O,  traveltimeasc) 

75 

■»6  *  Run  the  Calculator  of  Total  Time 

77  import  subprocess 

70  subprocess .call ( ’C;\\Work_Folder\\Clar)c_UH_Tot alTime_FAC.exe ' ) 


21,5.1,  Con(L0_0  “  22,10.2,  Con(LU_0 
42,1.6,  Con(Ln_0  =»  43,1.6,  Con(LU_0  = 
82,2.0,  Con(LU  0  =  90,1.6,  Con(LD  0  = 


=  23,15.3,  Con(LU_0  “  24,20.4, 
52,7.0,  Con(LU_0  =  71,7.0, 
95,1.6,  2.0)))))))))))))) 


*0efault  Workspace  (have  to  create  this) 
t  Overwrite  pre-existing  files 

IThis  is  sisqply  a  folder  for  the  output  data 
*outlet  location  of  the  basin 

IDEM  of  the  Area  (larger  than  the  actual  basin  size) 

IFAC  of  the  Area  (larger  than  the  actual  basin  size) 

IFlow  Dir  of  the  Area  (larger  than  the  actual  basin  size) 
*Land  Cover  of  the  Area  (larger  than  the  actual  basin  size) 
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Appendix  B:  C++  Code  to  Compute  TAH 


1  #include  <stdio.h> 

2  #include  <sCring.h> 

3  tinclude  <ntath.h> 

4  void  main ( void) t 

5  inc  r,  c,  ro,  co,  rol,  col,  ro2,  co2,  ncols,  nrows,  max_fac; 

6  double  h,  distance,  val,  xll,  yll,  cellsize,  ND,  t=60.0; 

^  char  JONKtl024]; 

8  bool  qpirt=false; 


9 

10 

11 

12 

13 

l-J 

15 

16 
17 
1£ 
19 
2Ci 
21 
22 
23 


//File  Names 

char  LN_Name [10241,  OF_Name [  10241 ,  rr_Name[10241,  Folder_Name{l='’C;\\Wor)c_Fclder\\"; 
strcpy  (LN_Name,  Folder_Naine) ;  strcat  (LN_Name,  "fac_a3c.asc") ; 
strcpy(OF  Name,  Folder  Name);  strcat (OF  Name, "totaltime . asc") ; 
strcpy  (TT__Name,  Folder_Neime) ;  strcat  (TT_Name,  "traveltaice  .asc") ; 

FILE  *LN,  *OF,  *rr; 


//Open  Input  and  Output  Files 
if(  (LN=fopen(LN_Name,  "r"))  !=  NULL  )[ 
elset  printf  (*'\n\nCould  Not  Open  %5,\n 
if(  (TT=fopen(TT_Name,  "r”))  !=  NULL  )t 

else!  printf ("\n\aCould  Not  Open  %s,\n 
OF=fopen(OF_Name,  "w") ; 


printf ("\n\nOpened  %s",  LN_Name) ; 

Please  Start  O'/ER!’!",  LN_Name  ) 
printf ( "\nOpened  %s",  TT_Name) ;  ) 
Please  Start  OVER!!!",  rT_Name  ) 


scanf ( "%s", 
scanf ("%s". 


&JUNK> ;  return;  ) 
&JUKK) ;  return;  ) 


24  //Read  and  Write  Header 


25 

fscanf (LN, 

"%s" 

"%d". 

£JUNK,  £ncols) , 

;  fscanf (TT, 

.  "%s"  "%d". 

£JUNK, 

£ncols) 

;  fprintf (OF 

,  "%s 

%d\n". 

JUNK,  ncols) 

26 

fscanf (LN, 

"%s" 

"%d", 

£  JUNK, 

£nrows) . 

;  fscanf (TT, 

"%5"  "%d". 

£  JUNK, 

£nrow5 ) 

;  fprintf (OF 

,  "%3 

%d\n". 

JUNK,  nrows) 

2~ 

fscanf (LN, 

"%s" 

"%lf", 

.  £JUNK, 

£xll) ; 

fscanf (TT, 

"%s"  "%lf’'. 

£  JUNK, 

£xll) ; 

fprintf (OF, 

"%s 

%lf\n", 

,  JUNK,  xll); 

2S 

fscanf (LN, 

"%s" 

"%lf", 

,  £JUNK, 

£yll); 

fscanf (TT, 

"%s"  "%lf". 

£JUNK, 

tyii); 

fprintf (OF, 

"%s 

%lf\n", 

,  JUNK,  yll) ; 

29  fscanf(LN,  "%s"  "%lf",  &JUNK,  Scellsize);  f3canf(TT,  "%s"  "%lf",  &JUNK,  Scellsize) ;  fprintf(OF,  "%s  %lf\n",  JUNK, 

30  fscanf(LN,  "%s"  "%lf",  &JUNK,  &ND) ;  fscanf(TT,  "%s"  "%lf",  &JUNK,  &ND) ;  fprincf(OF,  "%s  %lf",  JUNK,  ND) ; 

31 

32  //Allocate  memory  and  read  xn  the  data 

33  printf ("\nReading  through  the  Raster  FilesVn") ; 

34  float  **T  *  new  float* [nrows] ;  //Total  Time  (min) 

35  double  *’*R  =  new  double*  [nrows]  ;  //Rate  of  Time  (ztin/m) 

36  int  **D  =  new  inc* [nrows];  //Flow  Accumulation  Value 

37  for(ro=0;  ro<nrow3;  ro++) [ 

35  T[ro}  =  new  float [ncols] ;  R[ro)  =  new  double [ncols] ;  D[ro]  =  new  lnc[ncols]; 

39  fQr(co=0;  co<ncol3;  co++)  1 

40  T[rol (co]=-1.0; 

41  fscanf (TT, "%lf",&R[rol [co] ) ;  f scanf (LN, "%lf", &val) ;  D[ro] [co]  *  inc (val+0. 01) ;  )  > 

42 


cellsize) ; 


43  //Go  through  the  Eaoh  Cell 

44  printf  ("NnLcoJcing  ac  each  ceil  and  analyzing  the  Total  TimeNn"); 

4d  loriro— u;  rc<nrows;  ro++) i 

45  iprincf  (OF,  'Nn") ; 

47  for(co=C;  co<ncDls;  co-+) 1 

43  if  (C'(ro]  [col  <0)  (  T[rol  (co  1=0 . 0;  fprintf  (OF,  "%f  ",-£999.0);  continue;  ) 

43  quit=false;  val=0.0;  discance=0.0;  rol=rc;  col=co; 

sn  wh1  le  ( rpii  r.=f  alse)  • 

51  //Look  Diagonals  First 


55 

57 


maK_£ac  ~  C{rol] [col] ; 

r=(rol);  c=(col);  if(r>=0  &&  c>=0  &&  c<ncols  &&  r<nrow3Mif(  D(rl  (cj  >^iiax_fac;  {max_fac=D{r]  [c] ;  ro2=r;  co2=c;}} 
r=(rol-l);  c=(col-l);  iC(x>=0  aa  c>=C  c<ncols  aa  r<nrows){ir(  Dir) [c] >=max_rac) <max_rac=D(r J [c] ;  ro2=r;  co2=c; } • 

r=(iol-l);  c=(col+l);  if(r>=0  £&  c>=C  £&  c<ncols  &&  r<nrows){if(  D[r J [c] >=max_fac) {max_fac=D[r j [c] ;  ro2=r;  co2=c;)' 

r=(rQl+l);  c=(col-l);  if(r>=0  £&  c>=C  &&  c<ncol2  &&  r<nrow3){if(  D[rl [c] >=max_fac) (max_fac=D[rl [c] ;  ro2=r;  co2=c;}' 

r={rol+l);  c=(col+l);  if(r>=0  &&  c>=C  &&  c<ncols  £&  r<nrow3)[if(  D(rl[c]>=max  fac) (max  fac=D[rl [c] ;  ro2=r;  co2=c;}’ 


//Now  Look  at  Cardinal  Directions  (preferred  if  there  is  a  tie  between  this  and  a  diagonal  direction 


60 

r=(rol-l);  c=col; 

if (r>=0 

££ 

0=0 

££ 

c<ncol3 

££ 

r<nrowsi (if ( 

D(rl [c]>=max 

fac; (max 

_fac=Dcr]  [c] ; 

ro2=r; 

co2=c; ) ) 

€1 

r=rol;  3=(col— 1); 

if (r>=0 

e>=0 

cc 

c<ncol8 

r<nxowsi  (if ( 

Djr]  [c]  >=Riax 

fac'  {max 

fac»Drr) (3j ; 

ro2»r; 

co2=c; > } 

62 

r=rol;  c=(col+l) ; 

if (r>=0 

££ 

0=0 

££ 

ckncols 

££ 

r<nxows(  (if ( 

Djrl (cl>=max 

fac;  (max 

'fac=D[rj  [=1; 

ro2=r; 

co2=c; } ) 

63 

64 

1— (iul+1);  c— col; 

if 

&& 

O-O 

&& 

c<uculs 

££ 

x<xiiuws>  (iX  ( 

D|iJ  (c]y — luax 

Xac!  (oiax 

_Xac-D[iJ  [cl; 

iu2— i; 

coZ-c; ) J 

65  //Determine  the  distance  across  the  cell 

65  if(  rol— ro2  |i  col=co2  )(  h=1.0;  )  //This  means  the  adjacent  cell  is  in  a  cardinal  direction 

67  else(  h=sqrt(2.0);  }  //This  means  the  adjacent  cell  is  at  a  diagonal 

6J  h  =  h*csllsize;  distance  =  distance  +  h; 

63 

70  if(  colr»=ncols  ||  rol>=iirows  1)  coKC  ||  roKO  )  {quit=true.*} 

71  //This  IS  the  source  cell,  so  this  is  the  end  of  the  li.ne 

■^2  else  if  (D[ro2 ]  (co2| “~E(roll  [col]  )  f val  —  val  +  R(rol]  [col]  *h;  quit“-rue;  > 

73  //This  cell  has  already  been,  analyzec,  so  just  use  it's  value 

else  ir(  T[ro23ico2l  >  o.o  )(  vai  *  vai  +  0 .5*n* iR(roij (coij +Riro2| {co2 J )  +  T[ro2j[cc2i;  quit  =  true;  } 
75  else]  val  =  val  +  0 . 5*h* (R[rol) (col}4R[ro2 1 (co21 ) ;  col=co2;  rol=ro2;  }  } 

75  r[rcl[co)  =  val;  fprincf (OF, "4f  ",T[rol (co) ) ;  }  ) 

77  fcloseiLN);  fclosepF);  fclose(TT); 

7?  return;  ) 


NOTE:  The  contents  of  this  technical  note  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an  official 
endorsement  or  approval  of  the  use  of  such  products. 
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